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ABS TRACT 


The equation of motion for an elastically mounted, !inearly damped 
circular cylinder subjected to a time-dependent transverse fluid force 
has been formulated. The steady-state and transient solutions have 
been obtained in terms of the in-phase and out-of-phase force coeffi- 
cients, amplitude ratio, and the interna! damping of the cylinder. The 
results have shown that a cylinder may be se!f-excited in the range of 
Strouha!l frequencies close to the natural frequency of the cylinder. 

It has further been shown that the heuristic models previously used are 
inadequate to correctly predict the in-phase or the inertia component 
of the exciting force. Additional experiments and analysis could en- 


able one to predict the behavior of cables subjected to ocean currents. 
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NOMENCLATURE 


amplitude of cylinder oscillation 

drag coefficient as defined by equation (29) 
inertia coefficient as defined by equation (30) 
Fourler-averaged drag coefficient 
Fourier-averaged inertia coefficient 
diameter of the cylinder 

instantaneous total force acting on the cylinder 
period of oscillation 

Time 

instantaneous velocity 

maximum velocity in a cycle (U. = 27 A/T) 
mean velocity of the uniform flow 
dimensionless amplitude ratio 

dimensionless frequency parameter 
dimensfonless period parameter 

Reynolds number (Re = VD/%) 

Strouhal number (S = f0/V) 

vortex-shedding frequency 

cylinder oscillation frequency 

phase angle 

fluid kinematic viscosity 

fluid density 

cylinder density 


density ratio, P./D. 
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» INTRODUCTION 


Experience has shown that elastic structures of one or more degrees 
Of freedom, such as tal! slender strucTures, cables, etc., can extract 
energy from the flow about them and can develop catastrophic flow-induced 
oscillations. The understanding of this energy-extraction process is of 
paramount importance if one is either to eliminate or minimize if or To 
design the elastic structure such that it can withstand the osci|/ ations 
under the contemplated environmental! conditions. 

EVident!y, the real origin of the self-excited oscillations lies tn 
the basic mechanism of the vortex-shedding process. Consequently, the 
prediction of such oscillations requires a tractable theoretical-formula- 
Tion of the vortex-shedding phenomenon. The fact that such a formulation 
does not exist even for separated flow about a stationary cy!inder shows 
the complexity of the treatment of flow about oscillating bodies. The 
problem would have been relatively simpler had the vortex shedding and 
the forces arising from it been independent of the motion of the body. 
Experiments show that vortex generated forces give rise to oscillations, 
and the oscillations modify the forces, and so on, finally resulting in 
a steady-state oscillation such that the Strouhal frequency of the vortex 
shedding nearly coincides with the natural frequency of oscillation of 
the body. 

Mathematical models of flow-induced vibrations of bluff bodies and 
the response of circular cylinders to vortex shedding have been aptly 
described by Parkinson LI], and Currie, et al. [2] and will not be 


repeated here. 





It appears that among the various models considered so far, the 
"wake-oscillator'™ model of Hartlen and Currie [3] attracted considerable 
attention among the students of vibration analysis. This model predicts 
some of the observed features of the sel f-excited oscillations even 
though the reasons as fo why an oscillator of the van der Pol Type 
should at all describe the self-excited oscillations of an elastic body 
are not yet clear. There are several coefficients in the model which 
must be determined experimentally. There are, at present, wide varia- 
tions in the methods of measurement of some of the coefficients. For ex- 
ample, the damping of the vibrations of a given beam or cylinder is 
determined in still air or water depending on whether the self-excited 
oscillations take place in air or water. The added mass is assumed To 
be equal to the displaced mass. of the cylinder. There is no distinction 
between the material damping as would be determined in vacuum, and the 
damping caused by the fluid. 

It became clear, during the course of the present investigation, 
that the forces acting on a cylinder oscillating harmonically transverse 
TO a uniform flow should be measured and the in=phase and out-of-phase 
components of the force should be determined. This phase of the inves- 
tigation has been carried by Meyers L4] and Fortik [5] under the 
direction of Professor T. Sarpkaya, and the drag and inertia coeffi- 
cients for the lift force have been determined through the use of 
Morison's [6] equation. 

The purpose of the work reported herein is to make use of the 
measured force by combining it with the equation of motion of an 


elastically mounted, linearly damped, circular cylinder and to attempt 


10 





to predict the self-excited oscillations of the cylinder for given and 
directly measurable physical parameters such as the mass, length, 
diameter, material damping, and the natural frequency of the body, and 
the density and velocity of the fluid. 

lt will be assumed that the measurements and the predictions are 
in the same Reynolds number range and that the force coefficients do 
not measurably depend on the Reynolds number in the range under considera- 


meON « 


a 





Piemonte to lS 


A. FORCING FUNCTION 
The transverse force acting on a circular cylinder undergoing har- 
monic oscillations in the transverse direction to a uniform flow is 


expressed in terms of the Morison's equation as L4] 


UT 
CL = aa = oy mn? (—) (eis sin at 
=z 9, LDV VT 
UT 
canoes) Neos SU itl cos: t+ C1) 
oe = T T 


where eI and Cay are the inertia and drag coefficients, respectively. 
Evidently, the normalized force as well as the drag and inertia co- 
efficients depend on the parameters U_T/D ety Dee andr D/ Vin Neal ison 
parameter is directly proportional to the relative amplitude, and the 
second parameter may be regarded as the frequency parameter (D/VT) or 
the reduced velocity (VT/D). 

For the purpose under consideration it is more advantageous to ex- 
pand the |cos wt|cos wt term in series and to retain only ieent lesa 
term. This procedure yields, 

aC 2c (H)* sin wt - & cy? cos wt (2) 

VT VT 
Equation (2) could be simplified further by combining it with the equa- 


tion of motion for an elastically mounted cylinder. 


Ee 





B. ELASTICALLY MOUNTED CYLINDER 


Fig. 1. Elastically mounted, linearly 
damped cylinder. 


The equation of motion for the elastically mounted, linearly damped, 


and periodically forced cylinder may be written as (see Fig. |) 


me + cx + kx = F = Divan ic (3) 


| 
De L 
where m represents the mass of the cylinder, c the linear damping, and 


k the spring constant. The derivatives of x are taken with respect to 
Z 





Feal time Tt, i.e., o% s x, etc. 
dt 
Defining, 
er Oiemmee= 2 40°. tT =u t , x/D = x_., 
m n n ‘es 
ee a ney = 
C= > a? we eis SV/# DO ee A/D 
es Ming f ox 
Q= f/f, ato, xX = 
a 2m w in ; : a 


pas 





equation (3) may be reduced To 


ee 2 ee Z Pu QS Zz e 
X + 20%, - a aw TC TT (20K) Sin OT 
8 ZS eZ 
a gC, (20x) a” cos QT | (4) 


Defining p. = P./0. and noting that CoS panacea equation (4) re- 


duces To 


F : om Z, k l6 
Xx + 20x a Kes px 8 (Cy Sin QT ow Xe Cay EOS WT) Gey, 


or by a change of variable from a TO a x/A, and further simplifica- 


tion one has 


We t 2by e ves 0 MC Simic = — X C d | COST) (6) 
ei 

Equation (6) may be written as 

F. : a ee 

ee + ZY. 5 es Se sie cvr —' co) (7) 





in which 
= = P (8) 
and 
6A Ca 
acre ten ——>—__— Sot 
iw Co 


The coefficient EQ* may also be expressed in terms of the maximum 


lift coefficient defined by 


—‘) ; 
P, VIDE or Sees — CL) (10) 


Nj 
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Inserting equation (10) {nto equation (3) and simplifying one has 


| Z 


26 VeDLic) 
59¢ oe Cine 
2 
m Aw 
Nn 
S) ws 
(9° ope 
mer 53 ~~ (A/D) a 
27 


Combining equations (12) and (8) and simplifying one has 


(C S 





) 
Lom 27 Z Hous ae Z 
= C + (— 5 =C_),) } 
(A/D) VT.2 m| 3q77 ral) 


(>? 


[+ should be noted that if (A/D) Cy were independent of A/D, then 


| 
it is easy to show that (C) )/ (A/D) would also be independent of A/D 
Since the experiments have shown that Ca is independent of A/D. This 
would make the results independent of A/D, a conclusion which is contrary 
To observations. This proves that (A/D) Ca and (C) ),/6A/D) are indeed 
dependent on A/D in @ complex manner. But this dependence is rather 


weak, That its why slight changes in, for instance, Pp. can lead to large 


emenges in Cc, or %. 


C. SOLUTION OF THE EQUATION OF MOTION 


As shown above, a convenient form of the equation of motion is 
pees 2Gy + y.= Ent Sin (St = a) (14) 
‘9 r r 


|. Transtent Solution 


The Homogeneous equation corresponding to equation (14) is 


ee 2CY ey = 0 (15) 


is 





[ts characteristic equation reduces To 


Eero (16) 


whose roots are given by 


. Se ye Cy) 


Considering only the underdamped case, it follows that t < Il. 
Therefore, according to Boyce and Di Prima [7], the solution to the 
corresponding homogeneous equation, that is, the transient solution is 


given by 


Cen er sin (yl - c2 «t= 8) (18) 


eh nD 


where Yh and 8 will be determined from the [Atria condi | ons. 


ee steady-State Sela 1on 


The particular solution to equation (14) is a steady-state 
oscillation of the same frequency of the excitation. Thus, as indicated 


in Thomson [8], the particular solution is assumed to be of the form 


er = 5 sin L(Qt - a) - 6] . (19) 


where @ is The phase of the displacement with respect to the excitation, 


as shown in Fig. 2, and is given by 


— eee ee we 





Fig. 2. Vector diagram. 
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and 


arc ae (20) 


=oa 


o 
MN 


and is is given by 


2 
en (21) 


Poe Jc y-n2)2 + (2e0)? 


Adding the transient and the steady-state solutions, the following com- 


plete solution is obtained 


a - @ St sin (vl-t* t+B) + uy Sin (QT-a-¢) (22) 


Also, its derivative with respect to Tt Is given by 


Y= aecon iG Meee se 2 Vloc2 14a) t 


r 


+ us 2 cos (QT-a-o) (72) 


Saeyying the first initial condition that t = 0 corresponds To es Oo. 
one finds 
Y sin (-a-¢) (24) 
ee ee 2 

h sin B 
Applying the second initial condition that t = 0 corresponds to Me = 0, 


One obtains 


CES Sea /\-c* cos g) + 9 cos(-a-$) (25) 


Inserting the expression for Y_ given by equation (24) in equation (25), 


h 


iy, 





one has 


Ji? 


a cor(=a=o)- au 


Ben walge shan o( 
Also, recalling that as T >, x > -A sin 8, and that Ves x/A > -sin 6, 
it is seen from equation (22) that ‘5 =-|. Introducing this value in 


equation (21) and noting that E > 0, one obtains 


J(1-22)2 + (272)2 


92 


ae (27) 


Substituting the expressions for if and yy Into equation (22), one ob- 


Tains the following complete solution of the equation of motion 


= ane ee sin (vl-c? t+8) -sin (Qt-a-9) (28) 


The results obtained through the use of equation (28) will be 


presented in connection with the discussion of results. 
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tte ooUoo1ON Or RESULTS 


The time dependent force acting on the cylinder was expressed, as 


indicated previously, in terms of Ca and C Assuming this force to 


d|° 


be applied to an elastically mounted and linearly damped cylinder, the 
equation of motion for the cylinder was developed and solved in terms 
of the given. parameters and experimentally obtained coefficients. Two 


of these coefficients, C_, and C inertia and drag coefficients respec- 


m| ee 


tively, were introduced in the expression of the Time dependent force 
acting on the cylinder. Using the experimental apparatus and procedure 
described in L4], data were obtained, which together with the results 


from [4], were used to plot Cay and C,, as a function of VT/D for 


fyie= 0.25, 0.50, 0.75 and 0.84, as shown in Figs. 3 and 4. 
The results plotted in Fig. 3 show that for the region of interest, 
meee tor V1/D in the range of 4 to 7, and for the precision of the 


experiments, C_, is independent of A/D. 


m | 


Qn the other hand, C is dependent on both VT/D and A/D as shown 


d| 


in Fig. 4. As noted previously, (A/D) C,, and (C,) A/D) are not con- 


d| 


stants. Furthermore, it can be shown that (C) ) / CA/D) decreases as 


A/D increases. Figures 5, 6, 7 and 8 show (C as a function of VT/D 


Lm 
momen/D = 0.25, 0.50, 0.75 and 0.84 respectively. 

Therefore, an iterative procedure is needed to determine the complete 
solution of the equation of motion for a specific case, as indicated 
below. 


m@instly, for a given VT/D, C. lseeordinedai on mid. 5. hen, for 


a specific relative amp!itude A/D and the same VT/D, Cay is obtained from 


ag 





Fig. 4, and the value of (C Feseobmcainea trom Figs. 5,6, 7 or 8, oF 


er 
alternatively from equation (13). Then, E is calculated using equa- 
tion (12), since O.. is known. Substituting This Uaie neh E in equa- 
mom (27), a value for G is calculated, since Q = (aly (ED A plot 
of |1/E versus 2 for specific values of ¢ is shown in Fig. 9. Then, 

the values for a, and 8 are calculated from equations (9), (20) and 
(26), respectively. Substituting these results into the equation of 
motion and plotting the results, one obtains curves similar to those 
shown in Figs. 10, || and 12. These figures show clearly the effect of 
C on the response of the system. A lowly damped system results in an 
initial overshoot of the amplitude before it reaches a steady-state. 

A highly damped system, on the other hand, reaches its steady-state 
gradually from its initial rest position. 

Evidently, one can perform various parametric studies on the etfect 
of damping, density ratio, etc., on the transient and steady-states of 
the oscillation provided that sufficient experimental data are avail- 
able to obtain accurately enough the drag ne ene coefficients. IT 
is also evident that the practical problem is the determination of 
whether a cylinder of given natural frequency, density ratio, internal 
damping, etc., will be self excited in a fluid of given velocity. 

The determination of this problem requires an iteration not only 
over VT/D but also over A/D in order to obtain the appropriate values 
of the drag and inertia coefficients. Preliminary attempts to resolve 
this aspect of the problem have resulted in the computer program given 
in Appendix A. This program needs to be improved one additional data 
particularly to the lower and higher values of A/D where the force 


transfer coefficients vary most rapidly with A/D. 


Za) 





Then, it will be possible to improve the program to predict with 
greater precision the range of self-excited oscillations as well as 
their amp!litude. 

A computer program for the complete solution given by equation (28) 
and for the plotting of the results is given in Appendix 8. 

It is evident from The foregoing discussion that the prediction of 
the characteristics of the self-excited oscillations is based on the 
measured time-dependent transverse force rather Than on an assumed 
variation of the lift coefficient through The use of a heuristic model 
[3]. It was, therefore, desirable to compare the components of the 
alternating force obtained from The experimental data with those pre- 
dicted by the Hartlen and Currie L3] model. For this purpose the in- 
Phase and out-of-phase components of the Hartlen and Currie transverse 
force coefficients have been written equal to the corresponding com- 


ponents of the mode! used herein Lsee equation (2)]. This procedure 


yielded 
S27 (a) 
C (Hartlen and Currie) = ~ —————C (29) 
d ~ d | 
Lye 
D 
and 
20° (5) 
C_(Hartlen and Currie) = a ar | (30) 
Viele a 


(—) 


in which Gee and Cay are the drag and inertia coefficients evaluated in 
The present investigation. 


Equations (29) and (30) have been evaluated for A/D = 0.25 and 0.50 


and the results are plotted as a function of oe er bse | osandal 4. 


Ze 





A similar plot obtained from [3] is shown in Fig. 15 for A/D = 0.2. 

Even though the A/D values are not exactly identical, Figs. 13, 14 and 
15 could certainly be compared as far as the overal|! variation of Fe 

and C are concerned. It is evident that C curves are fairly similar 
and that the Hartlen and Currie model could be adjusted sufficiently to 
yield similar C values for identical A/D values. However, Cc. values 
which are comparable only for ee < | show opposite variations. It 
has been previously stated that the use of an added mass coefficient 

in terms of the displaced mass of the body is not sufficient and that 
the inertia coefficient of a body oscillating in the transverse direc- 
tion to a uniform flow could be significantly different from its dis- 
placed mass. It suffices to conclude that the Hartlen and Currie mode! 
does not predict Ce in accordance with the observations and that the 
overall predictions of the model may be attributed partly to the adjust- 
ment of the other free parameters such as the damping ratio, the density 


Ponto, etc. 
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IV. CONCLUS IONS 


The results presented herein have shown that (a) a cylinder can 
undergo self-excited oscillations for certain values of the ratio of 
the natural frequency of the cylinder to the Strouhal frequency of the 
vortices (f/f), internal damping of the cylinder (ft), and the density 
ratio (0); (b) the drag and inertia coefficients must be determined 
experimentally as a function of A/D and VT/D and must be incorporated 
into the equation of motion; (c) the prediction of the sel f-excited 
oscillations and their amplitude will be possible with the acquisition 
of additional data for the transverse force; and that (d) the results 
presented herein are quite encouraging and certainly a reflection of 
the combination of the fluid dynamical and vibrational aspects of the 


phenomenon unlike the heuristic models previously proposed. 
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Normalized amp|itude of the JIft force versus VI/D for A/D 
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Normalized amplitude of the 11 ft force versus VT/D for A/D = 0.50. 
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Normalized amp|itude of the lift force versus VT/D for A/D = 0.84. 
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